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HAMILTONIAN  DYNAMICS  OF  A 
RIGID  BODY  IN  A  CENTRAL  GRAVITATIONAL  FIELD 


Li-Sheng  Wang*  P.S.  Krishnaprasad*  J.H.  Maddocksf 

ABSTRACT.  This  paper  concerns  the  dynamics  of  a  rigid  body  moving  under  the  influence  of  a 
central  gravitational  field.  Explicit  account  is  taken  of  effects  arising  because  of  the  finite  extent 
of  the  body.  The  hamiltonian  framework  of  the  problem  is  exploited  to  elucidate  questions 
concerning  approximation,  symmetry,  Poisson  reduction,  relative  equilibria,  and  associated 
stability  problems. 


1.  INTRODUCTION. 

In  the  study  of  the  Newtonian  (gravitational)  many-body  problem,  it  is  customary  to  treat 
the  bodies  as  point  masses.  See  (Sternberg  [40],  Smale  [39],  and  Abraham  and  Marsden  [1]). 
However  the  proper  accounting  of  stable  planetary  spins  for  instance,  would  seem  to  require  the 
consideration  of  rigid  (possibly  nonhomogeneous)  bodies  of  finite  extent  as  a  first  approximation. 
The  works  of  Duboshin  [7],  Ermenko  [10],  Elipe,  Ferrer  and  Cid  [8],  [9],  are  concerned  with 
the  existence  of  special  solutions  (e.g.  central  configurations)  in  the  Newtonian  many-rigid-body 
problem.  However,  in  these  papers,  the  natural  geometric  and  group-theoretic  underpinnings  of 
the  problem  are  not  exploited  to  the  extent  possible.  We  are  not  aware  of  a  systematic  program 
along  these  lines. 
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In  the  design  of  large  earth  satellites,  aerospace  engineers  have  had  to  account  for  a  “gravity- 
gradient  torque”  and  its  effect  on  the  stability  of  earth-pointing  satellite  attitude.  In  the  literature 
related  to  this  problem,  there  are  studies  of  relative  equilibria  and  quasi-periodic  motions  based 
on  various  approximate  models  of  the  coupling  between  orbital  motion  and  attitude  motion  of 
earth  satellites.  We  refer  the  reader  to  the  work  of  Beletskii  [2],  Duboshin  [6],  Roberson  [34], 
[33],  Longman  [22],  Meirovitch  [27],  Mohan,  Breakwell  and  Lange  [28],  Likins  [21],  Sincarsin 
and  Hughes  [38],  Pascal  [30],  and  Sarychev  [36].  The  basic  problem  at  hand  is  the  dynamics  of 
a  rigid  body  or  gyrostat  in  a  central  gravitational  field. 

In  the  present  paper,  we  work  out  the  noncanonical  hamiltonian  structure  of  the  problem 
of  motion  of  a  rigid  body  in  a  central  gravitational  field.  The  group  , S' 0(3)  of  three  dimensional 
rotations  appears  as  a  symmetry  group.  Poisson  reduction  by  the  action  of  50(3)  yields  a  nine¬ 
dimensional  system  that  corresponds  to  observing  the  dynamics  from  a  moving  frame.  In  this 
body  frame,  the  dynamics  manifests  the  effect  of  a  fictitious  torque  known  as  the  gravity-  gradient 
torque.  There  are  Casimir  functions  that  are  conserved  independently  of  the  hamiltonian  and 
hence  of  any  (convenient)  approximations  to  the  Newtonian  potential.  We  shall  compute  relative 
equilibria  and  determine  their  stability.  All  motions  (whether  exact  or  approximate)  remain 
confined  to  the  level  sets  of  Casimir  functions  which  are  eight  or  six  dimensional  symplectic 
leaves. 

It  is  noteworthy  that  the  Poisson  structure  for  the  finite  dimensional  problem  studied  here  is 
closely  related  to  the  one  used  by  Krishnaprasad  and  Marsden  [17]  in  their  study  of  the  dynamics 
of  a  rigid  body  with  a  flexible  attachment  (a  physically  distinct  and  infinite-dimensional  problem). 
A  key  link  is  the  geometry  of  Poisson  reduction. 

2.  CONFIGURATION  SPACE. 

In  Figure  1,  let  C  denote  a  fixed  gravitating  body  of  mass  M  (with  spherical  symmetry) 
that  influences  the  motion  of  a  rigid  body  B  of  mass  m .  The  inertial  frame  of  reference  (of 
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the  observer)  is  attached  to  C  and  a  body  frame  is  fixed  on  the  rigid  body  13  at  its  center 
of  mass.  A  typical  material  particle  Q  in  the  rigid  body  is  represented  by  the  inertial  vector 
q  =  DQ  +  r,  where  B  is  an  element  of  50(3)  (independent  of  the  particle)  and  r  is  the 
vector  from  C  to  the  center  of  mass  of  body  B.  At  any  instant,  the  configuration  of  the  rigid 
body  B  is  determined  uniquely  from  the  pair  (B,r)  e  SE( 3),  the  special  Euclidean  group  of 
rigid  motions  in  It3  . 

In  what  follows,  we  will  see  that  this  is  an  example  of  a  simple  mechanical  system  with 
symmetry  in  the  sense  of  Smale  [39]  (see  also  Abraham  and  Marsden  [1]).  Appendix  1  includes 
a  short  introduction  to  the  abstract  framework. 


Figure  1.  Rigid  Body  in  Central  Force  Field 


3.  LAGRANGIAN. 

The  kinetic  energy  of  the  rigid  body  relative  to  the  observer  at  C  is, 
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T  =  \  J  I  ?  I2  dm(Q) 

B 

where  dm  (•)  denotes  the  mass  measure  of  the  body.  Here  onwards,  |  •  |  denotes  the  Euclidean 
norm  in  It3 .  It  is  an  elementary  fact  that  the  above  expression  simplifies  to  the  formula 

T  =  \  <0,m>  +  j  |  r  |2 

where  f l  is  the  body  angular  velocity  vector  of  the  rigid  body  ,  m 
body  and  I  is  the  moment  of  inertia  tensor  of  B  in  the  body  frame 

Recall  that  the  body  angular  velocity  f l  is  defined  by 

B  =  BCl, 

where 

0  — 

0  =  0  — 

—  ^1  0 

is  the  skew-symmetric  matrix  associated  to  fi . 

The  spatial  angular  velocity  u  is  defined  by 

B  =  uB, 

and  we  have  the  relation 

c o  —  BSl. 

In  the  notation  of  Appendix  1,  we  note  that  K 
SE( 3)  the  configuration  space. 

The  gravitational  potential  energy  of  the  body  B 


=  2 T  defines  a  riemanman  metric  on 


is  given  by, 


(3.1) 

is  the  total  mass  of  the 
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(3.2) 


V=  -  J 

B 

f  GM  dm(Q ) 

J  | r  +  BQ | 


where  G  is  the  universal  gravitational  constant. 

The  Lagrangian  for  the  problem  is  then  a  function 


L  :  T  (SE( 3))  -+  1R, 

(P,r,Q,r)  T  -  V". 


4.  SYMMETRY. 

The  inertial  observer  at  C  has  the  freedom  to  change  his  frame  of  reference  to  a  new 
orientation.  This  corresponds  to  an  50(3)  action  on  the  configuration  space  C  =  SE(3 )  : 

$:  50(3)  xChC 

(P,(R,r))  (PB,Pr). 

It  is  easily  checked  that  this  action  leaves  invariant  the  kinetic  energy  T  (  riemannian  metric 
on  C  )  and  the  potential  V. 

The  hamiltonian  H  =  (T  +  V)  is  given  by, 

H  =  I  <n,  I-1  n>  +  hh  _  /  (Q),  (4.D 

B 

where  II  =  I  is  the  body  angular  momentum  of  the  rigid  body  B ,  and  p  =  mr  is  the 
spatial  linear  momentum  of  the  body.  One  has  also  the  formula, 
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7T  =  Bin  +  r  X  mr 
=  BII  +  r  X  p 

for  the  spatial  angular  momentum  of  the  rigid  body. 

It  can  be  verified  that  tt  —  tc  (n,B,r,p)  is  an  Ad*-equivariant  momentum  mapping  for 
the  lifted  action  on  T* SE( 3)  and  hence  is  a  conserved  quantity  for  the  dynamics  Xjj  . 

This  is  further  equivalent  to  Euler’s  balance  law.  To  see  this,  let  Resultant  denote  the  force 
resultant  on  the  rigid  body.  Then, 


P 


X resultant. 


(4.2) 


On  the  other  hand,  the  torque  resultant, 


^  f  (r  +  BQ)  x  (r  +  BQ)  GM  , 

•^resultant  —  —  I  - j - ; „  „  ,, - dm\Q) 

J  |  r  +  BQ  ]3 

B 

=  0. 

Thus  angular  momentum  (or  Euler’s)  balance  law  yields: 
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Collecting  together  the  balance  laws  one  can  write  the  spatial  form  of  the  dynamics  as 


Equivalently,  in  mixed  body  and  space  variables  (II,  5,  r, p)  we  get: 

GM(BTr  x  Q) 


n  =  n  x  r 


w 

B 


r  +  BQ 


dm(Q), 


P 


=  ~J 


B 


GM  (r  +  BQ) 
I  r  +  BQ  |3 


dm(Q), 


B=  BI-'Il, 
r  =  p/m. 


(4.5) 


Now,  since  H  is  50(3) -invariant,  one  can  induce  a  hamiltonian  H  on  the  quotient 
T*  (SE( 3))  /50(3)  and  express  the  dynamics  X ^  in  terms  of  appropriate  reduced  variables 
(see  Appendix  2  for  the  general  framework).  In  the  present  context  it  is  easy  to  determine  the 
reduced  variables.  Note  that 

:  50(3)  x  T*SE{ 3)  ->  T*SE( 3) 

(R,(U,B,r,p))  (II,  RB,  Rr,  Rp) 

is  the  lifted  action  on  T*SE( 3).  A  representative  for  each  equivalence  class  in  T* SE(3)/ SE(3) 
is  given  by 


(n,  /,  BTr,BTp). 
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Thus  the  reduced  variables  (or  convected  variables)  are: 

II,  the  body  angular  momentum,  A  =  BTr ,  the  convected  radius  vector  from  C,  and 
fi  =  BTp,  the  convected  linear  momentum. 

In  terms  of  these  convected  variables,  the  dynamics  Xfj  takes  the  form 


and  the  hamiltonian  H  is  given  by, 

»  =  +  j  (4'7) 

B 

Equations  (4.6)  with  hamiltonian  (4.7)  are  the  Poisson  reduced  equations  on  T*  S  E(3)  /  SO(3) 
~  so* (3)  X  IR3  X  IR3  .  In  terms  of  the  Poisson  tensor  A  on  so*  (3)  X  1R3  x  1R3  derived  in  Appendix 
2,  these  equations  take  the  compact  form, 

/li\  ffl  A  A  /Vn^\ 

A  =  A  0  I  VXH  =  A  VH.  (4.8) 

V A/  \ A  -I  0 J  \VpHj 

The  Poisson  structure  is  rank-degenerate,  and  there  are  nontrivial  Casimir  functions  of  II,  A,  ji . 

(  Casimir  functions  are  kinematic  conserved  quantities  for  equations  of  the  form  (4.8).  See 
Appendix  2  for  the  precise  definition.  )  In  fact,  any  function  C $  of  the  form 

CV  =  <H|n  +  Ax/r  I2), 

is  a  Casimir  function.  Here  <j)  :  IR  — »  IR  is  any  smooth  scalar  function.  Moreover,  these  are 
the  only  Casimir  functions  defined  on  the  open  set  of  generic  points  of  A  . 
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From  the  general  properties  of  Casimir  functions  (see  Appendix  2)  we  know  that  is 

an  integral  invariant  for  any  hamiltonian  vector  field  and  in  particular  for  .  It  is  further 

important  to  note  that  replacing  H  by  a  suitable  approximation  (such  as  derived  from  series 
expansions  of  the  Newtonian  potential  term)  does  not  affect  the  integral  invariance  of  C This 
is  of  some  use  in  developing  an  analytic  perturbation  theory. 

5.  RELATIVE  EQUILIBRIA. 

The  concept  of  relative  equilibrium  goes  back  to  Poincare.  For  simple  mechanical  systems 
with  symmetry,  there  is  an  elegant  characterization  of  relative  equilibria  due  to  Smale  [39].  We 
discuss  this  below  and  use  it  in  computing  relative  equilibria  for  a  rigid  body  in  a  central  force 
field. 

Let  (M,  A",  V,  G )  be  a  simple  mechanical  system  with  symmetry  as  defined  in  Appendix  1. 
Assume  that  the  lifted  action  of  G  acts  on  T*M  freely  and  properly.  Then  the  quotient  space 
T*M/G  is  a  smooth  manifold  with  an  induced  Poisson  structure.  Let  {•,  -}o  be  the  canonical 
Poisson  bracket  on  T*M.  Given  /,  g  G  C°°  (T*M/G),  the  induced  Poisson  bracket  of  / 
and  g  is  defined  by 


{/,  9}  0  f  =  {f°f,goT}0 

where  r  :  T*M  — »  T*M/G  is  the  canonical  projection.  For  any  G -invariant  hamiltonian 
function  H  on  T* M ,  we  have  the  induced  function  H  :  T*M/G  — >  1R  defined  by, 

H  o  f(x)  —  H(x). 

In  terms  of  the  induced  Poisson  structure,  and  H ,  the  projected  hamiltonian  vector  field  Xfj 
on  T*M/G  is  defined  by  the  condition,  for  any  /  e  C°°(T*M/G), 
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XM  =  {/,  H}. 


DEFINITION. 

ze  e  T*M  is  relative  equilibrium  for  H  if 

('f(Ze))  =  0. 

For  the  dynamics  of  a  rigid  body  in  a  central  gravitational  field,  the  relative  equilibria 
are  determined  by  setting  the  time  derivatives  in  equation  (4.6)  (or  (4.8))  to  zero.  On  the 
other  hand,  in  general  position,  i.e.  II  ^  0,  VC^  spans  the  kernel  of  A.  Thus  we  have  the 
energy- Casimir  characterization  of  relative  equilibria  in  general  position:  (II,  A,  ff)  is  a  relative 
equilibrium  iff 

XH  =  VC*,  for  suitable  <j>,  (EC) 

iff  (  Lagrange  multiplier  characterization  ) 

i-xn\  /  n  +  Axfi  \ 

VaF  =  c  /i  x  (n  +  A  x  ff)  .  (LM) 

l  )  \(II  +  Ax^)xA/ 

where  c  ^  0  is  a  constant  and 

^ (A)  =  -  X 

Before  we  proceed  to  solve  (EC),  we  note  an  alternate  characterization  of  relative  equilibria. 
REMARK  5.1. 

It  can  be  shown  [l]  that  ze  is  a  relative  equilibrium  iff  there  exists  a  (  6  3  such  that  the 
flow  of  Xh  , 

FXu  (ze)  =  exp(t£)  (*e), 

(i.e.  the  dynamical  orbit  is  simply  a  group  orbit).  Thus  if  the  observer  were  to  be  set  in  uniform 
motion  according  to  the  one-parameter  group  exp(tff) ,  then  for  such  a  moving  observer,  a  relative 
equilibrium  will  appear  to  be  stationary. 
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In  the  work  of  Smale  [39],  there  appears  a  characterization  of  relative  equilibria  for  simple 
mechanical  systems  with  symmetry  (see  Appendix  1  for  the  relevant  notation).  We  present  this 
as  an  algorithm  to  determine  relative  equilibria. 

Algorithm  for  relative  equilibria 

Step  0.  Pick  £  e  Si,  the  Lie  algebra  of  G  .  Let  £ m  denote  the  corresponding  vector  field 
on  M  determined  by  the  action  $  of  G . 

Step  1.  Search  for  the  critical  points  qe  e  M  of  the  function  (  the  augmented  potential  ) 

Vt  :  M  -+1R 

V(M  =  V(q)  -  ~  K  «„(<?),  £„(«))  (5.1) 

Step  2.  Find  the  corresponding  conjugate  momentum  pe  by  the  formula 

Pe  =  Ii'  ««(«.))•  (5.2) 

Then  =  (qe,  Pe)  is  a  relative  equilibrium  for  the  hamiltonian  function 

H(c,)  =  r(9)  +  i/c((A-‘)-1c„(A'‘)-1a5). 

This  principle  was  used  in  [41]  to  determine  relative  equilibria  for  the  dynamics  of  two  rigid  bodies 
connected  by  a  ball-in-socket  joint.  In  what  follows  we  apply  this  principle  to  find  the  relative 
equilibria  for  the  problem  of  rigid  body  motion  in  a  central  force  field.  Here,  the  elements  of  a 
simple  mechanical  system  with  symmetry  consist  of 
M  =  SE(  3), 

K((U!,  U2),  (Wi,  W2))  =  tr(U\  I'  W[r)  +  m<U2,  W2  >e, 

V(B’  T)  =  ~  L  |  r  TbQ  I  dm(Q)’ 

and 

G  =  50(3), 
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where  <  •,  •  >e  denotes  the  Euclidean  inner  product  on  1R3  , 

(Uu  U2),  (Wu  W2)  6  T(B,r)SE{ 3), 

and  I '  is  the  coefficient  of  inertia  of  the  rigid  body.  The  superscript  T  in  Wi  denotes  matrix 
transpose. 

For  £  e  so(3),  the  corresponding  infinitesimal  generator  of  the  group  action  on  M  can  be 
found  as 

r)  =  ( £b ,  |r). 

We  then  have 


K(Cm(B,  r),  £m(B,  r))  =  <  BT£,  I BT  t  >E  +  m\£  x  r 


and 


Vf(B,  r)  =  -  i  <  Bt(,  I  BTi  >B  -  1  m  U  X  r  |2 

”  L  I  r  f  %Q  |  dmm  (5'3) 

We  then  get  the  first  order  conditions  for  ( B ,  r)  to  be  a  critical  point: 

(0  m  Z  x  (f  x  r)  +  ^  =  0 

(«)  t x  (bibt£)  -  dm(Q)  =  o. 

Next,  we  calculate  pe  in  Step  2.  The  map  K ^  can  be  found  as  follows.  For  (wiB,  W2) , 

(v,iB,U2)  €  T(B,r)SE{ 3), 
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Thus 


A'b  (wiB,  w2)  ( u\B ,  u2 ) 

=  tr(wiB  l'BTu^)  +  m  <  w2 ,  u2  >e 

rp 

—  <  u i,  AIA  uq  +  <  u2,  mw2  >e  ■ 


K\wxB,  w2)  =  ((_ BIBTWl )  B,  mw2)  G  T(B  r)  SE( 3). 


We  then  have 


pe  =  A'b(f„(,«))  =  ((BIBT£)  B,  m|r) 

=  (B(lS?0,  r)  . 

Note  that  in  the  forimila  for  pe  ,  the  two  components  correspond  to  the  angular  momentum 
and  linear  momentum  respectively.  If  we  let  ft  denote  the  body  representation  of  the  linear 
momentum,  we  get 


(Hi)  ft  =  mBT£r. 

Substituting  f l  =  B1  £,  A  =  BTr,  conditions  («),  (n),  (in)  read 

(A  +  Q)GM 


(i')  m  ft  x  (ft  x  A)  +  j 

Jb 

(n)  nxm  -  f 

Jb 


B  \  ^  +  Q  \3 
(A  x  Q)  GM 


dm(Q)  =  0, 


b  I  ^  +  Q  !3 


dm(Q)  =  0, 


and 


(in')  ft  =  m(ftx  A). 


(5.4) 


These  conditions  are  identical  to  the  conditions  obtained  from  the  reduced  dynamics  (4.6)  and 
the  definition  of  relative  equilibrium. 
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Now,  if  we  take  the  cross  product  with  A  on  both  sides  of  (*'),  we  get 


m\e  x  (fic  x  (Qe  x  Ae))  +  /  X  Q^  Gfdm(Q)  =  0. 

Jb  I  A  e  +  Q  r 

(Here  again  the  subscripts  e  refer  to  equilibrium.)  Comparing  it  with  (u'),  we  obtain 

m  Ae  x  (fle  x  (fle  x  Ae))  +  x  I^e  =  0. 

By  standard  identities  in  vector  analysis,  we  get 

fle  x  (I  —  m\e\T)  Qe  =  0.  (5-5) 

We  conclude  that  £le  must  be  an  eigenvector  of  the  matrix  I  —  mAeA^. 

Let  ke  denote  the  corresponding  eigenvalue.  Then  one  can  obtain  the  relative  equilibrium 
characterization  (LM)  from  (5.4)  by  setting, 

1 

ke  +  m  |  Ae  |2 

Conversely,  using  the  identity, 

A  x  ((fi  x  A)  x  O)  =  —  H  x  (A  x  (fi  x  A)) 

=  (A  •  Sl)n, 

and  a  few  further  algebraic  manipulations,  one  can  derive  (5.4)  from  the  relative  equilibrium 
characterization  (LM).  We  leave  the  verification  to  the  reader.  Thus  the  two  characterizations 
are  equivalent.  Of  course,  for  simple  mechanical  systems  with  symmetry,  the  equivalence  of  the 
energy-Casimir  characterization  (EC)  or  (LM)  and  the  variational  characterization  based  on  the 
augmented  potential  holds  in  general. 
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Note  that  we  fix  £  while  searching  for  critical  points  of  .  Thus  ft  =  Br £  is  of  fixed 
norm  as  B  varies  over  SO (3). 

Let 

I  «  |2  =  u  l2  =  /?• 

Define  Vp  (ft,  A)  to  be  V^(B,  7)  expressed  in  the  convected  variables  ft,  A.  Then, 

~  1  777 

Vp  (ft,  A)  =  -  -  <  ft,  m  >E  -  -~\a  x  a  |2 

+  V  (A).  (5.7) 

Clearly,  the  critical  points  of  Vp  on  the  sphere  |  ft  J2  =  f3  satisfy  the  unconstrained  variational 
principle, 

We  +  ~  I  «  I2 )  =  0.  (5.8) 

where  1/2 c  is  a  Lagrange  multiplier.  The  first  order  conditions  associated  to  (5.8)  are, 

I  ft  +  m  A  x  (ft  x  A)  =  —  ft  (5.9a) 

m  (ft  x  A)  x  ft  =  VA  V.  (5.9 b) 

These  are  exactly  the  equations  we  get  by  eliminating  p  =  m(ft  x  A)  in  the  relative  equilibrium 
characterization  (LM).  The  unconstrained  variational  principle  (5.8),  parametrized  by  c,  and  the 
associated  first  order  conditions  (5.9)  appear  to  be  most  suited  to  the  explicit  computation  of 
relative  equilibria.  Before  we  proceed  with  such  specific  computations  we  make  some  general 
geometric  observations  concerning  relative  equilibria. 
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Observe  that,  by  taking  the  inner  product  of  both  sides  of  {in')  with  Ae ,  we  get 


<  Ae,  fle  >E  —  0  (5.10) 

at  a  relative  equilibrium  (  Ae,  fie,  )•  If  (  re,  Bt  )  is  a  relative  equilibrium  configuration,  then 
the  dynamical  motion  is  such  that 

r(t)  =  e**  re 

B(t)  =  e**  Bt.  (5.11) 

This  follows  from  Remark  (5.1)  that  at  a  relative  equilibrium  the  dynamical  orbit  is  just  a  group 
orbit. 

PROPOSITION  5.1. 

In  relative  equilibrium,  the  radius  vector  r(f)  generates  a  right  circular  cone. 

Proof 

From  (5.11), 

<  r{t),  r(t )  >e  =  <  re,  re  >E  =  <  re,  re  >E  ■ 

Also 


< 

<  r,£  >£ 

£  <r,Z>E  z 

l£l2 

lei2 

=  < 

<r,Z>E 

r,r>E-  |f|2 

=  < 

1 

63 

A 

s- 

<  re,  £  >\ 

l£l2 

=  < 

re,  re  >E  - 

<  re ,  £  >% 

lei2 

=  constant. 

.  2  .  i 

Thus  r(f)  is  a  circle  of  radius  (<  re,re  >E  —  <rj^p>g  J  centered  at  C'  =  <r\tfE  £  •  See  Figure  2. 

| 
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NON-GREAT  CIRCLE  MOTIONS. 

For  a  rigid  body  of  finite  extent,  if  the  center  of  (relative  equilibrium)  rotation  C'  does  not 
coincide  with  the  center  C  of  the  force  field,  then  the  stationary  motion  will  be  called  a  non-great 
circle  motion.  The  existence  of  such  motions  is  in  question.  See,  e.g.  the  model  problem  below 
and  also  the  gyrostat  example  in  Rumyantsev  [35]. 


From  equation  (4.5), 

GM  (r  +  BQ) 

B 


-/ 


r  +  BQ  |3 


dm(Q )  —  p  =  m  r 


d 2 


=  m- 


■it*  (e‘'r0  =  me''?'2r" 
Substituting  r  =  et(-  re  and  B  =  Be  on  the  left  hand  side,  we  get, 


f  GM  (re  +  Be  ^dm(C)  -  m  G  r 

J  1  re  +  BeQ  |3  dm{Q) 


Taking  the  inner  product  of  both  sides  with  £,  we  get, 


/ 


GM 


<  C  re  >E  Be  Q  >E 

I  re  +  BeQ  I3 


dm(Q )  =  m  <  £,  re  >  ==  0. 
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Hence 


<  £,  re  >E 


f  dm(Q)  f 

<  f,  BeQ  >E 

J  re  +  BeQ  |3  ~  J 

I  re  +  BeQ  |3 

dm(Q). 


B 


B 


Equivalently, 


<  re  >e 


-s 


<£,  Be  Q  >e 
|  re  +  BeQ  |3 


dm(Q)  /  J 


dm(Q) 


re  +  BeQ  |3 ' 


The  quantity  <  £,  re  >e  is  proportional  to  the  cos(#)  (refer  to  Figure  2),  and  C  and  C'  coincide 
iff  <  £,  re  >E  =0.  If  the  body  B  were  a  point  mass,  Q  =  0  and  hence  <  £,  re  >e  =  0.  If  for 
a  rigid  body  of  finite  extent,  the  integral 

<  Be  Q  >E 


1 


re  +  BeQ  |3 


dm{Q)  ^  0, 


then  C,  C'  are  not  coincident. 
Since 


<  re  >E  =  <  Bjt,  Bj  re  >e 
—  f2e  •  Ag, 

we  conclude  that  a  relative  equilibrium  (Ae,  Q,e,  /ie)  determines  a  non-great  circle  motion  iff 


Oe  •  Ae  ^  0.  (5.12) 

One  can  test  the  non-vanishing  condition  (5.12)  in  various  settings.  We  now  demonstrate  that 
there  are  examples  which  do  not  admit  great  circle  relative  equilibria.  We  first  assume  that  the 
relative  equilibrium  is  a  great  circle.  Then  the  equilibrium  can  be  found  by  solving  (  from  (5.9) 
and  $2  •  A  =  0  ), 

si  x  m  =  o, 

/»W*"W)  =  m|afi  (5'13) 
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with  m|fl|2  being  the  Lagrange  multiplier. 

Now  we  consider  a  model  problem.  The  body  is  an  asymmetric  “molecule”  consisting  of  six 
point  masses,  two  on  each  principal  axis.  See  Figure  3. 
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In  this  example,  we  know  that  I  is  diagonal  and  thus  for  a  great  circle  solution,  must  he 
along  one  axis.  For  a  given  set  of  data,  for  example, 

mxi  =  101,  mX2  =  1,  my±  -  100,  my2  =  1,  mz\  —  99,  mz 2  —  1, 

=  0.01,  Vl  =  0.01,  zx  =  0.01, 

the  corresponding  x2,  y2,  z2  can  then  be  determined  such  that  (0,  0,  0)  is  the  center  of  mass.  The 
corresponding  function  V"t  on  the  sphere  |A|  =  400  is  shown  in  Figure  4.  The  coordinate  system 
is  the  following.  The  sphere  is  parametrized  by  the  usual  spherical  coordinates  9,  0  <  9  <  ir ,  and 
<j),  0  <  <f>  <  2tt.  In  Figure  4,  the  function  is  plotted  above  the  disc  of  radius  7 r,  with  (9,<f>) 
interpreted  as  planar  polar  coordinates. 


Figure  4.  The  function  V t 


The  extremal  critical  points  are  determined  numerically  to  be  as  follows. 

maxima  :  A  =  (-398.5,-33.7,-7.2),  (399.3,-22.1,-10.6), 
minima  :  A  =  (13.7,32.4,398.5),  (-4.0,8.5,-399.9). 


20 


These  extremals  are  found  by  using  an  optimization  package  CONSOLE  [11]  with  the  assistance 
of  the  3-D  graphical  representation  in  Figure  4.  From  Figure  4,  we  note  that  there  are  also  two 
saddles  near  the  axis  t<i .  We  therefore  search  for  the  maxima  along  that  axis  and  check  if  they 
are  also  minima  in  the  transverse  direction.  By  this  process,  we  can  verify  that  no  critical  value 
of  A  is  perpendicular  to  a  principal  axis.  Accordingly,  we  conclude  that  for  this  example,  there 
are  no  great  circle  relative  equilibria. 


6.  APPROXIMATIONS. 

For  typical  applications  in  the  modeling  of  planets  or  artificial  earth  satellites,  the  nominal 
radius  of  the  orbital  motion  is  very  large  compared  to  the  dimensions  of  the  orbiting  body. 
Accordingly,  it  seems  appropriate  to  consider  various  approximations  of  the  gravitational  potential 
based  on  Taylor  series  in  a  neighborhood  of  |  A  |  =  oo  or  equivalently  |  r  |  =  oo.  Whilst  such 
approximations  are  common  in  the  literature  cited  in  the  Introduction,  it  is  unclear  whether  the 
symmetries  and  conservation  laws  inherent  in  the  problem  are  respected  by  the  approximation 
process. 

In  the  present  paper,  we  take  the  Poisson  reduced  model  (4.8)  as  the  logical  starting  point 
for  approximations.  The  hamiltonian  H  is  approximated  to  various  orders  of  e  =  (nominal 
dimension  of  body)  /  (orbital  radius),  by  the  Taylor  series  expansion  of  the  V  (A)  potential  term 
appearing  in  H : 


f  GM 

V{X)  =  -LiTTQ 

GM 


dm(Q ) 


^  I  JB 

GMm 

TTT 


dm(Q )  jl  — 
GM 


<  Q,  A  > 
I  A  I2 


1  1  Q  l2 

2  I  A  P 


3  <Q,  A  >s 


+  o(\  AT4) 


+ 


2|  A 


,tN  3  GM  x  x 

tV  ^  +  2  |Tp  <  A’  1  A  > 


2  |  A  |4 

+  o(i  a  n 


(6.1) 


In  (6.1)  the  first  term  in  brackets  is  of  the  order  e°  and  the  next  term  is  of  the  order  e2  .  The  e1 
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term  is  absent  due  to  the  vanishing  of  fB  Q  dm(Q). 

We  will  therefore  consider  two  approximate  model  Hamiltonians, 

Ho 

and 

H2 

Upon  substituting  Hq  and  H2  respectively  for  H  in  the  Poisson  reduced  dynamics  (4.8),  one 
obtains  the  order  zero  reduced  dynamics; 


and  the  order  two  reduced  dynamics; 


As  already  noted  at  the  end  of  section  4,  all  such  approximations  admit  a  common  set  of  conserved 
quantities  (Casimir  Functions)  of  the  form  C<$,  =  <j>  (|  II  +  A  X  /i  |2) .  Since  the  order  0  dynamics 


i<n,  i-n>  +M--  ^ 


-<n,i-‘n>  + 


GMm 


A  ±  ,T,  3  GM 

|3  tr ®  +  2  I A  I5 


<  A,  IA  > 
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is  essentially  decoupled,  it  has  additional  conserved  quantities  of  the  form  ip (|  II  |2) ,  and  the  spin 
energy  1/2  <11,  I— 1 II  > .  If  the  body  is  spherically  symmetric,  i.e.,  I  =  kl ,  then  the  order  two 
approximation  collapses  to  the  order  zero  approximation.  In  general,  the  order  two  approximation 
displays  nontrivial  spin-orbit  coupling. 


6.1  ORDER  ZERO  RELATIVE  EQUILIBRIA. 

With  the  order  zero  approximation  of  H ,  the  relative  equilibria  (fle,  \e )  satisfy 

Jc  fig 

m(f2e  x  Ae)  x  fle. 


I  fie 
GMm 


a  F  Ae 


(6.6  a) 

(6.6  b) 


By  taking  the  inner  product  of  both  sides  of  (6.6b)  with  fle ,  we  conclude  that  Ae  •  fle  =  0,  i.e. 
all  relative  equilibria  in  the  order  zero  approximation  give  rise  to  great-circle  orbits.  From  6.6(b) 
and  the  condition  fle  ■  Ae  =  0 ,  we  get  the  Kepler  frequency  formula, 

(  GM  \1/2 


fl  I  = 


(6.7) 


Summarizing,  the  only  relative  equilibria  for  the  order  zero  approximation  are 


(a)  fle  is  a  principal  axis  of  I; 

(b)  Ae  is  a  vector  perpendicular  to  fle  satisfying  the  Kepler  formula  (6.7);  (6.8) 

(c)  fxe  =  m  (fle  x  Ae)  completes  a  triad. 


With  the  same  assumptions  as  (6.8a)  and  (6.8b)  above,  it  is  possible  verify  the  existence  of 
“uniformly  spinning  solutions”  to  the  order  zero  reduced  dynamics: 


n(0  = 
m  = 

/F)  = 


m(t)=  me 


U) 


exp  (  t  fie  )  Ae 


OJ  \  I  LO 

m  (  1  +  r  I  exP  ( r  )  fie  Ae, 


fie 


fie 


(6.9) 


23 


with  the  modified  Kepler  frequency  formula, 


(u  -fi  |  |) 


(sr 


(6.10) 


The  quantity  uj  measures  the  body  spin  relative  to  a  moving  Frenet-  Serret  frame  at  the  center 
of  mass  of  the  body. 

6.2  ORDER  TWO  RELATIVE  EQUILIBRIA. 

The  first  order  conditions  for  the  variational  principle  (5.8)  take  the  form 


(I  —  m  AAr)0  =  (-  —  m  |  A  |2)  f2 

c 


(6.11) 


m  |  |2  A  —  m  (Q  •  A)  = 


GMm 


3GM 
21  A  p 


<r(I)A 


3 GM  15  GM  .  T  ... 

+  m? IA  -  v  mr  1  A>A 


The  equations  (6.11)  admit  a  family  of  solutions  (relative  equilibria)  corresponding  to  great  circle 
motions: 

(a)  Qe  is  a  principal  axis  (eigen-vector)  of  I  with  corresponding  principal  moment  of  inertia 

i  =  1,  2,  3; 

(b)  Ae  is  a  principal  axis  (eigen- vector)  of  I  perpendicular  to  Vte ,  with  associated  principal 
moment  of  inertia  J;-; 

(c)  /j,e  -  m  (Qe  x  Ae); 

and,  the  following  modified  Kepler  frequency  formula  holds: 
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In  the  above  relation  *,  j,  k  are  distinct  and  takes  values  in  {1,  2,  3}.  Hence  the  correction 
term  in  (6.12)  may  be  of  either  sign.  It  follows  that  for  the  order  two  approximation  there  are 
twenty-four  1  parameter  families  of  relative  equilibria  (accounting  for  12  being  in  each  of  the  six 
directions  parallel  to  the  principal  axes  (with  sign)  and  four  directions  for  A  corresponding  to 
each  choice  of  12),  the  scalar  parameter  being  /?  —  |  12  |2  =  |  £  |2  as  in  Section  5. 

This  conclusion  appears  to  be  a  classical  result  exhibited  in  different  form.  See  for  instance 
the  book  of  Beletskii  [2].  However,  the  hamiltonian  point  of  view  together  with  the  approach  of 
reduction  has  entirely  eliminated  the  formidable  mess  of  Euler  angles  and  such. 

In  the  following,  we  show  that  for  practical  parameter  ranges,  all  the  relative  equilibria  in 
the  order  two  approximate  model  are  great  circle  motions.  Let 

(I  —  mAAr)12  —  a:f2, 
or 

I  12  -  o!2  —  mAAr12. 

With  the  notation  r  =  A  ^12,  we  have 


mr  A  =  112  —  a  12.  (6.13) 

We  note  that  r  ^  0  corresponds  to  solutions  that  are  not  great  circles,  while  r  =  0  implies  a 
standard  eigenvalue  problem.  The  dot  product  of  (6.13)  with  12  then  yields 


M 


-(12TI12  —  mr 2), 


o 


and  substitution  in  (6.13)  gives 


-  ^pfinP-oa^ia  + 


(6.14) 


Taking  the  dot  product  with  0,  of  the  second  equation  in  (6.11),  we  get  the  following  equation 


3_ 

2|AP 


1 5  o 

(m  +  ^trl~  2P^aTia)t  +  TViTfiTlA  =  °- 


|A| 


(6.15) 


Assuming  r  ^  0,  and  multiplying  (6.15)  by  r,  we  have 


(m  +  ^trl 


15  \Tt\\_2 


2|A|4AJIA)T“  +  |A|2^TItA  0- 


With  the  expression  for  rA  in  (6.14),  we  obtain  the  equality, 


3  15  .T--.  3 


(m+2iIp'rI'  2|A|4 


ATlA+ |A|2|fipflTlS!^2  m|A|2|S2|2 


{|0|2|lfl|2  -  |flTIO|2}  .  (6.16) 


But  we  know  that 


|0|2|io|2  -  |oTm|2  >  o, 


thus  (6.16)  can  have  a  solution  with  r  ^  0  only  if 


q  15  3 

m  +  2\)^ trI  ~  2|A]4  ATlA  +  |A|2|0|2 ^ ^  ~  °’ 


which  can  be  true  only  if 
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m 


15 

2|A|4 


AJ 


IA  <  0, 


or 


15 

2m 


FIX  >  l, 


(6.17) 


where 


A  =  ,  I  =  J  j||ydm(Q)  1  -  J  ~^din{Q). 

It  is  easy  to  see  that  for  large  A ,  (6.17)  is  not  satisfied.  In  particular,  the  ratio  must  be  greater 
than  Ys  •  But  for  typical  artificial  satellites,  this  ratio  is  approximately  10~10  .  For  the  motion  of 
moon  around  the  earth,  it  is  approximately  1.6  x  10~5  .  Thus  we  have  shown  that  for  the  practical 
case  of  large  orbit  radii,  the  24  relative  equilibria  (for  the  second  order  approximate  model )  are 
the  only  relative  equilibria.  This  conclusion  is  of  special  interest  since  we  have  constructed  a 
numerical  example  (  the  “molecule”  in  Section  5  )  in  which  the  exact  model  has  no  great  circle 
relative  equilibria. 


7.  STABILITY  OF  RELATIVE  EQUILIBRIA  IN  THE  APPROXIMATE  MODELS. 

In  this  section,  we  study  the  stability  properties  of  the  relative  equilibria  for  the  approximate 
models  discussed  in  Section  6.  For  both  cases,  the  triple  (n,  A,p.)  is  a  relative  equilibrium  if  the 
three  vectors  are  along  the  three  principal  axes.  Without  loss  of  generality  we  let 


ne  —  |n|  ei  —  Ii|0|ei, 

Ae  =  |  A  |  e2, 


l^e  — 


mlH||A| 

h 
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m|ft|A|  e3, 


(7.1) 
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where  \Q\  and  |A|  are  related  through  the  appropriate  form  of  the  Kepler  frequency  formula,  and 

I  ei  =  I\ei, 

I  ^2  =  he-2: 

I  e3  =  -^363  ■ 

We  shall  examine  the  stability  of  this  relative  equilibrium  in  various  cases  determined  by  the 
relative  magnitudes  of  the  principal  moments  of  inertia  . 


7.1  Order  Zero  Approximate  Model.  (  Instability  Proof  ) 

For  the  order  zero  reduced  dynamics  (6.4),  the  energy-Casimir  method  of  Appendix  3  is 
inconclusive  since  the  second  variation  of  the  energy-Casimir  function  is  only  positive  semi-definite 
(has  a  zero  eigenvalue).  We  linearize  the  system  around  the  relative  equilibrium  (7.1).  Let 

8x  =  (<5IIi,  8TI2,  <5113,  8X2,  8X3,  8fi  1,  8f.i 2,  8^)t. 


We  have  the  linearized  system 


8x 


A  8x, 


where 


/ 


A 
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0 
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By  the  frequency  formula  (6.7),  we  can  write  A  in  the  form 


( 

0 

0 
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Denote  the  upper  left  3x3  matrix  by  B  and  the  lower  right  6x6  matrix  by  C .  It  can  be  shown 
that 

Pi(s)  =  det(sl  -  B)  =  -s  (^.s2  +  llj  J1-2  , 

p2(s)  =  det(sl  —  C)  =  s2  (s2  +  |fl|2)2  • 

The  characteristic  polynomial  of  A  is  p(s)  =  pi(s)p2(s) .  It  can  be  further  verified  that  the 
minimal  polynomial  of  A  is 

m(s)  =  (^  +  l3~JlJ2~/lw2)  (a»  +  ima). 

The  occurrence  of  a  repeated  root  of  the  minimal  polynomial  at  s  —  0  implies  linear  instability 
of  the  relative  equilibrium  (7.1)  for  the  order  zero  approximate  model  (See  Gantmacher  [12], 
Theorem  3,  pp.  129).  Alternatively,  the  one  parameter  family  of  “uniformly  spinning  solutions” 
given  by  (6.9)  represents  a  perturbation  of  the  relative  equilibrium  (7.1)  that  departs  any  small 
neighborhood  of  the  relative  equilibrium  in  finite  time,  and  hence  we  have  instability.  We  note 
that  this  conclusion  is  independent  of  the  relative  magnitudes  of  the  I;’ s. 
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REMARK. 


The  projection  of  (II,  A,  yu )  to  the  space  of  II  projects  the  order  zero  dynamics  to  the  usual 
rigid  body  dynamics.  For  this  projected  dynamics ,  the  equilibria  in  which  the  vector  II  is  along 
the  maximum  or  minimum  principal  axes  are  stable. 

I 


7.2  Order  Two  Approximate  Model.  (  Energy- Casimir  Method  ) 

We  now  study  the  stability  of  relative  equilibria  of  the  order  two  reduced  dynamics  (6.5). 
For  the  relative  equilibrium  (7.1),  we  have  the  following  identity,  (  from  (6.11)  or  (6.12)  ) 


m|J7|2 


a  \  t 


m  + 


2|A| 


rtr(I) 


2|A|: 


(7.2) 


Now  we  discuss  sufficient  conditions  for  the  stability  of  these  24  relative  equilibria. 

In  general,  for  a  hamiltonian  system  with  hamiltonian  H  and  Casimir  C,  we  have  the 
energy-Casimir  type  sufficient  condition  for  stability  (see  Appendix  3).  The  general  form  of  the 
energy- Casimir  function  for  our  case  is 


Ha 


-nTi-1n  +  ; 


l|/r|2  mGM  GM  /T.  3 GMiTt,  ,/llTT  x  i2  . 

2  m  fXj  W'r(1)  +  WA  U  +  H5'n  +  '  X  1  ' 


The  first  variation  of  can  be  found  as  II,  A,yu)  =  \7H<f,  ■  6x:  where 

/  I^n  +  ^n  \ 


V^(II,  A,yu) 


GM  ( 

|A|3  v 


m  + 


3/rI  15A 1  IA 


V 


2 1 A  | 2  2 1 A  | 4 


A  +  ^jfpTA  +  fi'pAi 


jl  _  a'  : 

in 


and  Sx  is  as  in  Section  7.1,  and, 


n  =  II  +  A  x  yu. 


The  matrix  representation  of  the  second  variation  of  is 


30 


Ill  the  above  formulae,  <j)'  represents  its  value  at  |n|2/2,  and  the  same  convention  is  applied  to 
(f)" .  Now  we  find  the  variations  at  (ne,Ae,/xe).  By  using  (7.2),  we  have 

/  Mei  \ 

V^(ne,Ae,/xe)  =  (1  +  <f>'K)  m|A||^|2e2  , 

\  |A||fi|e3  / 

where 

K  =  h+m\X\2. 

Thus  in  order  for  the  first  variation  to  vanish,  we  require  1  +  (j>'K  —  0,  or 

<7-3) 

Substituting  these  values  in  the  second  variation  formula,  we  get  F  =  V277?)(IIe,  Xe,  jie)  as 
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where 


m|fi||A| 
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m|f2||A| 
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m|fJ  |  j  A  | 
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m |  Q 1 1 A  | 
'  R 


-m|0|2(4+^-)  0 

I  2  rn  G  M 
+  |A|> 


-mi(i  +  =^) 


fl-m|A|2 

mR 


-  -  t"K*\n\\ 


(7.4a) 


T-1  .  .*2  I\  3 GM  ,  ■, 

F44  =  m|S2|  —  +  ...5  (h  -  h), 


(7.46) 


T-i  i/-ai2  3 GM.  .  . 

Fee  =  m|0|  4 — jyp— (^3  —  h)- 


(7.4c) 


With  the  lower  triangular  matrix  L  defined  by 


m|f2||A|/i 

R-h 


/2m|n||A| 

K-I2 
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with 


Lgi  =  — (m|£7|£95  +  1)  , 

iL  —  1 1 

r  uM 

i9!  “  ft? 

0»  =  -mi  (i + f5^-) , 


D55 


m|flj2  (  4  + 


K~h\ 

R-h) 


2  mGM 

~FF 


we  can  transform  F  into  a  diagonal  matrix 


D  =  L  F  Lt, 


(7.5  a) 
(7.5  b) 
(7.5c) 
(7.5  d) 


where 


D  =  diag 


{ 


R  —  Ii 

hR 


K-I2 

hK  ’ 


A"  -  A 
/3a 


7  (7i  —  h)Di4c, 


—  D55-,  A66, 


h~h 

rn(K  -  I3)  ’ 


(J3  —  I2)D88, 


uh  R  -  A- 

D55  m(R  -  ii) 
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with 


D44 


m|0|2  3GM 

J7-72  +  "[aF 


7?88 


3GM 
m|A|5F66 ' 


(Since  congruence  transformations  preserve  the  matrix  inertia ,  we  can  read  off  the  number  of 
negative  eigenvalues  of  A  from  D .) 

We  shall  now  consider  the  case  in  which, 

I\  >  I3  >  7*2  • 
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To  have  stability  from  the  energy-Casimir  method,  we  require  that  all  the  eigenvalues  of  F 
(  equivalently  of  D  )  be  positive.  This  holds  if 


>  0, 


D  55  <  0, 


and 


For  (7.4a),  we  have 


V},  R-K 

D55  m(R  —  Ji) 


0. 


ra|A|2 

K 


+  <t>"I1K2\Cl\2. 


(7.6  a) 


(7-6  b) 


(7-6  c) 


Thus  (7.6a)  holds  if  <j>"  >  0.  We  next  consider  (7.6b).  By  the  definition  (7.5d)  of  D55  and  the 
frequency  formula  (7.2),  we  get 


D55 


K  -  h 
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1  +  2m]\\2't:*'I  2m^A|2^2} 


Let 
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2m|A|2 


(h  +  J3 


2/2). 


For  the  case  under  consideration  e  >  0  and  for  |A|  large  compared  to  the  typical  dimensions  of 

the  body,  e  is  small.  The  other  term, 

K-h  m\\\2  (1/K  -<f>"K2\Q\2) 

R-h  ~  m\\\2/K  +  (/>"hK2\tt\2  ' 


For  [A |  large  and  <f> "  large  enough, 

K-h  m\\\2  (-(/>" K2\Q\2)  m|A|2  a 

R-h  ~  <t>"hK2\Sl\2  h 


Thus 
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D5 5  ~  m|Oj2  (  4  - 


1  +  e 


e 


for  |A|  large  and  <f>"  large  enough.  Since  0  »  4,  we  have  D55  <  0.  This  is  (7.6b).  Now  we  look 
at  (7.6c).  It  is  easy  to  see  that  if  we  show 


tt2  R-  K  ^ 

^59  H - 7+ - TT-^55  <  0’ 

m(R  —  Ii) 


(7.8) 


then  together  with  (7.6b),  we  have  (7.6c).  From  (7.7),  for  |A|  large  and  <j>"  large  enough, 

U59  ~  -|ft|(l-(9). 


From  the  definition  (7.4a)  of  R,  we  have 

R-K 

R 


=  (f>"K3\n\2. 


Thus 
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<j>"K3 1  ft  |5 


m(R-h)  +  |Q|2  Ji 


A" 

-  —  =  1+6 , 


for  |A|  large  and  <f>"  large  enough.  Now  we  verify  (7.8).  Under  the  same  condition, 


vi,  +  *  im»(i  -  «)2  +  i(i  +  #)m|n|2  (4  2 


m 


1  +  e  / 


-  N2(=-T77+(1-rt7)e 

~  -|0|2^ 

<  0. 


Thus  (7.4)  hold  for  |A|  large  and  f"  large  enough.  We  have  the  following  theorem. 


STABILITY  THEOREM. 

For  the  order  two  approximate  model,  the  relative  equilibrium 
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in  =  Jin 


i  a  =  i2\ 

I  n  -  hn 


is  stable  if  |A|  is  sufficiently  large  and, 

I\  >  J3  >  12- 


I 

This  shows  that  the  relative  equilibrium  in  which  the  body  center  of  mass  traverses  a  circular 
orbit,  the  angular  velocity  lies  along  the  principal  axis  of  the  body  with  the  largest  associated 
moment  of  inertia  (  minor  axis  of  the  ellipsoid  of  inertia  ),  and  the  radius  vector  is  aligned  to  the 
principal  axis  with  the  least  associated  moment  of  inertia  (  major  axis  of  the  ellipsoid  of  inertia  ), 
is  a  stable  relative  equilibrium. 

HISTORICAL  REMARK. 

A  similar  theorem  appears  in  Beletskii’s  book  [2],  pp.  94-102.  Beletskii  uses  a  spatial/inertial 
model  of  the  coupling  between  translational  and  rotational  motion  and  presents  arguments  based 
on  a  Lyapunov-Chetayev  approach  [5],  and  uses  in  effect  the  variational  equations  about  the 
stationary  motion.  In  contrast,  here  we  make  consistent  use  of  modern  hamiltonian  methods  and 
reduced  variables.  The  methods  of  this  paper  yield  a  nonlinear  stability  theorem  and  generalize 
to  nonrigid  and  other  complex  configurations.  See  for  instance  the  examples  considered  in  [17], 
[23],  [16],  [18],  [19],  [32]- 

7.3  Order  Two  Approximate  Model.  (  A  Lagrange  Multiplier  Approach  ) 

The  previous  section  demonstrated  that  it  is  sometimes  not  straightforward  to  explicitly  find 
an  appropriate  function  <j)  in  the  energy- Casimir  method.  In  Appendix  3,  we  describe  a  more 
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classical  characterization  of  relative  equilibria  as  critical  points  of  the  constrained  variational 
principle, 


min  #2(11,  A,  /i)  (7.9) 

subject  to  C(Il,X,n)  —  constant 

where  H2  is  the  hamiltonian  (6.3)  and  C  is  the  Casimir  |jll  +  A  X  fi\2 .  The  associated  first- 
order  conditions  coincide  with  the  characterization  (LM)  of  relative  equilibria,  with  the  unknown 
constant  c  being  interpreted  as  a  Lagrange  multiplier. 

The  Lagrangian  ( in  the  sense  of  optimization  theory  )  associated  with  the  above  constrained 
variational  principle  is  recovered  if  in  Section  7.2  we  take  <j)(x)  =  —  c  x .  Consequently  the  second 
variation  can  be  recovered  as  a  special  case  of  that  calculated  in  Section  7.2.  When  (f>  is  linear, 
<j>"  =  0 ,  and  consequently,  R  =  I( .  Therefore  the  second  variation  of  H  —  c  C  reduces  to 
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A  comparatively  simple  Gaussian  elimination  then  reveals  that  Fc  is  congruent  to 
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Dc  =  diag 


K  -  h 
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where 


Dc 

-‘-''55 


5m|f2|2 


2mGM 


and  the  other  coefficients  are  as  defined  previously.  From  the  identity  (7.2),  we  have  further  that 

2  \ 


Dh  =  m\a\2  5  - 


1  +  € 


For  |  A  |  sufficiently  large,  the  expressions  for  the  various  coefficients  reveal  that  K  is  large 
and  positive,  and  -D44,  D£5  ,  Fqq  and  are  all  positive.  Consequently  the  signs  of  the  entries 
of  Dc  are  determined  by  the  signs  of  the  entries 


{  +  ,  +  ,  +  ,  (h  —  h)  ,  —  ,  +  ,  (h  —  h)  ,  (h  ~  h)  i  +  }  • 


We  shall  restrict  attention  to  satellites  in  which  the  inertias  are  distinct  so  that  Fc  is 
nonsingular.  Otherwise  additional  symmetries  arise,  and  the  analysis  is  slightly  more  complicated. 
There  are  six  cases  of  distinct  inertias.  See  the  following  table,  in  which  the  number  of  negative 
eigenvalues  of  Fc  is  shown  in  each  case,  and  each  case  is  assigned  a  reference  number  in 
parentheses. 


I 2  min 

I 2  middle 

I2  max 

>  h 

(1)  1 

(2)  2 

(3)  3 

h  <  h 

(4)  2 

(5)  3 

(6)  4 

Table  1. 

According  to  Theorem  A3. 2  (in  Appendix  3),  it  suffices  to  analyze  whether  the  condition 


(h,  Fch)  >0,  V  ft  #  0  and  (VC(IIe,  Ae,  ge),  h)  =  0 


(7.10) 
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is  satisfied.  Because  the  subspace  of  admissible  variations  h  has  codimension  1,  condition  (7.10) 
cannot  hold  whenever  Fc  has  two  or  more  negative  eigenvalues.  Accordingly  the  only  case  in 
which  (7.10)  might  hold  is  case  (1),  I\  >  1$  >  I2,  in  which  Fc  is  nonsingular  and  has  precisely 
one  negative  eigenvalue. 

To  analyze  (7.10)  in  case  (1),  we  shall  apply  a  general  result  bearing  upon  families  of 
extremals  to  variational  principles.  Notice  that  (7.1)  actually  defines  a  one-parameter  family  of 
relative  equilibria  which  can  be  regarded  as  being  parametrized  by  the  magnitude  of  the  radius 
of  the  orbit,  i.e.  |A|.  But  the  multiplier  c  is  related  to  |A|  through  c  =  ji+^n mr  (cf.  (7.3)),  so 
the  family  can  also  be  parametrized  by  the  multiplier  c.  Along  this  family  the  Csisimir  can  be 
written,  using  (7.1)  and  (7.2),  as 


■  |n  T  A  x  n I"* 


i|ft|2(/1+m|A|2)2 

GM(h  +  mjA|2)2 
m|A|5 


+  2^1  _  ^ 


(7.11) 


Consequently,  for  j  A  |  large  the  Casimir  is  an  increasing  function  of  ]  A  |  along  the  family  of  relative 
equilibria,  and  consequently  a  decreasing  function  of  c.  We  may  now  apply  the  aforementioned 
result. 


LEMMA.  (Maddocks  [23],  Lemma  5.2,  pp.  316) 

Suppose  a  family  of  variational  principles  of  the  type  (7.9)  have  a  family  of  critical  points 
xe(c)  parametrized  by  the  multiplier  c.  Moreover,  suppose  that  the  second  variation  at  a 
particular  extremal  is  nonsingular  with  one  negative  eigenvalue.  Then  the  second-order  sufficient 
conditions  (7.10)  at  that  extremal  are  satisfied  if  and  only  if  the  constraint  C  is  a  decreasing 
function  of  the  multiplier  c  at  that  parameter  value. 

COROLLARY. 

Solutions  (7.1)  in  case  (1),  h  >  h  >  h  are  Lyapunov  stable  for  all  |A|  sufficiently  large. 
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Proof 


It  has  been  shown  that  the  hypotheses  of  the  previous  Lemma  hold,  and  that  C(c )  is 
decreasing  for  relative  equilibria  with  |A|  sufficiently  large.  Thus  condition  (7.10)  holds  and 
Theorem  A3. 2  in  Appendix  3  then  applies  to  provide  the  desired  result. 

I 

Accordingly  we  have  rederived  the  Stability  Theorem  proved  in  Section  7.2.  As  a  final 
remark,  we  observe  that  we  have  not  proven  instability  in  the  order  two  model.  Actually  the 
results  of  Maddocks  [24] (Section  5)  can  be  applied  to  show  that  for  large  |A|  the  relative  equilibria 
in  any  of  the  cases  (2), (4), (6)  in  Table  1,  are  dynamically  unstable.  An  outline  of  the  analysis 
is  that  when  Fc  has  an  even  number  of  negative  eigenvalues  and  C(c)  is  a  decreasing  function, 
then  the  linearized  dynamics  must  possess  an  unstable  real  eigenvalue. 

8.  CONCLUSIONS 

This  paper  represents  a  first  step  in  our  program  to  understand  the  geometry  and  dynamics 
of  motion  in  a  central  gravitational  field.  Treating  the  rigid  body  as  a  model  problem,  we  adopt 
a  modern  approach  to  hamiltonian  mechanics  to  address  questions  concerning  approximation, 
symmetry  and  reduction,  Poisson  structure,  relative  equilibria,  and  associated  nonlinear  stability 
problems.  Our  methods  should  extend  naturally  to  problems  of  current  interest  such  as  the 
dynamics  of  tether-connected  (  and  other  )  multibody  systems  in  orbit,  the  dynamics  of  elastic 
shells,  etc.  A  common  thread  in  the  program  would  be  the  use  of  geometrically  exact  models 
(  i.e.  models  that  respect  natural  groups  of  symmetries  ). 

Among  the  new  results  of  this  paper  we  note  the  reduced  model(s),  examination  of  non¬ 
great  circle  solutions,  the  instability  proof  for  the  order  zero  approximation,  and  the  rigorous 
proof  of  a  nonlinear  stability  theorem  in  the  order  two  approximation. 

The  geometric  framework  adopted  here  should  be  helpful  in  understanding  and  further 
exploring  some  of  the  deep  and  exciting  questions  that  have  emerged  in  recent  years  such  as: 
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spin-orbit  coupling  and  the  problem  of  Hyperion  [43]  [13]  [14]  [31]  [35];  the  stable  resonances  of 
Markeev  and  Sokolskii  [25];  the  work  of  Beletskii  [3]. 

Further  questions  concerning  bifurcations  and  instability  are  of  some  interest  and  appear 
to  be  worthy  of  careful  study.  In  addition  to  the  methods  used  in  the  present  paper,  such 
investigations  may  be  facilitated  by  some  of  the  new  techniques  based  on  the  energy-momentum 
method  [37].  We  hope  to  report  on  these  problems  at  a  later  date. 
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APPENDIX  1.  (Simple  Mechanical  Systems  with  Symmetry) 

Let  (M,I<)  be  a  riemannian  manifold  and  let  G  be  a  Lie  group  with  associated  action, 

$  :  G  x  M  ->  M 

(g,q)  ^  ®g(v) 

where  in  an  isometry  for  each  g  e  G.  The  riemannian  metric  induces  a  vector  bundle 

isomorphism, 


Kh  :  TM  ->  T*M 


defined  by 


Kh  (vg)  ■  Wg  =  K  (vq,wq), 

for  all  vq:wq  e  TMq ■  The  canonical  symplectic  structure  uj  =  ~d90  on  T*M  can  be  pulled 
back  to 

n  -  (I&r  (w), 

which  is  also  an  exact  symplectic  structure  on  TM.  One  can  verify  that  the  action  lifts  to 
symplectic  actions  $T  and  on  TM  and  T*M  respectively. 

Let  V  :  M  —>  ] R  be  a  G—  invariant  (potential)  function  on  M.  The  hamiltonian 
H  :  T*M  — +  1R ,  is  defined  by, 

H  (a,)  =  i  A'((/C)-‘  a„  (/C)"1  a,  )  +  V  o  rj,  (a,) 
where  tm  ■  T*M  — >  M  is  the  canonical  projection. 


46 


The  hamiltonian  vector  field  Xfj  on  T*M  determined  by  H ,  is  given  by  the  condition, 


dH(Y)  =  w(Ih,  F), 


for  all  vector  fields  Y  on  T*M.  The  hamiltonian  system  ( T*M,lv,Xh )  or  equivalently  the 
quadruple  (M,  K,  V,  G )  is  a  simple  mechanical  system  with  symmetry.  It  admits  a  momentum 
mapping  in  a  natural  way.  To  see  this,  let  A  denote  the  Lie  algebra  of  G  and  A*  the  dual 
space  of  G.  The  symplectic  action  on  T*M ,  defines  a  Lie  algebra  homomorphism  of  A 

into  hamiltonian  vector  fields  on  T*M;  we  denote  this  correspondence  as  £  i— >  Then 

the  map, 

J  :  T*M  -»•  A* 


defined  by, 


J(<Xq)-Z  =  o)  (a9),  £eA 


is  an  Ad*  -equivariant  momentum  mapping.  Hence  J  is  a  conserved  quantity  of  the  system 
(T*M,lo,Xh)-  (See  Abraham  and  Marsden  [1]  for  proofs.) 
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APPENDIX  2.  (Poisson  Structure  and  Reduction) 

In  this  appendix  we  outline  the  essentials  of  Poisson  reduction  and  derive  the  Poisson  bracket 
applicable  to  the  problem  of  rigid  body  motion  in  a  central  force  field.  A  good  source  for  Poisson 
structures  is  the  book  by  Libermann  and  Marie  [20].  See  also  Olver  [29],  the  papers  of  Weinstein 
[42],  Marsden  and  Ratiu  [26],  and  Krishnaprasad  and  Marsden  [17]. 

A  Poisson  manifold  P  is  a  smooth  manifold  equipped  with  an  IR-bilinear  map  (Poisson 
structure)  on  the  space  of  smooth  functions, 

:  C°°(P)  x  C°°(P) C°°(P) 

satisfying  the  axioms,  for  /,  g  G  C°°(P), 

(i)  {f,g}p  =  -  { gJ}p 

(ii)  {fg,h}p  =  g{f,h}P  +  f{g,h}P 

(iii)  {/>  {g,h}p}P  +  {g,  {h,f}p}P  +  {h,  {/,  g}P}P  =  0. 

Associated  to  a  Poisson  structure,  there  is  a  unique  twice  contravariant  skew-symmetric, 
smooth  tensor  field  A  on  P  such  that 

{f,g}p  =  A  (df ,  dg), 

where  df ,  dg  are  differentials  of  / ,  g ,  respectively.  The  tensor  field  A  defines  a  vector-bundle 
morphism, 

A#  :  T*P  ->  TP 

ax  ^  A #{ax)  G  TXP 

satisfying, 

M A#K))  =  A  (*)(«,,/?*)  for  all  fdx  €  T*XP. 

Let  G  be  a  Lie  group  and  let  :  G  X  P  — >  P,  ( g ,  x  )  *-»•  ^s(^)>  be  a  group  action  such 

that  'kg(-)  is  a  Poisson  morphism  for  every  g  6  G.  Suppose  that  the  action  is  proper  and 
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free.  Then  there  exists  a  good  quotient  P/G  that  carries  a  Poisson  structure  {■, -}p/a  induced 
from  the  one  on  P  satisfying,  for  f,g  £  C°°(P/G) , 

{f,g}p/G°ir  =  {/0  7r,?0  7r}p. 

Here  tt  :  P  — *  P/ G  is  the  canonical  projection.  By  construction,  it  is  a  Poisson  morphism. 

(?-equivariant  dynamics  on  P  induce  dynamics  on  P/G.  Suppose  h  :  P  — >  H  is  a 
G -invariant  hamiltonian  function  on  P  ,  i.e., 

h(^g  (x))  =  h  (x)  V  g  e  G. 


Define  a  vector  field  Xh  by 


X„[f\  =  {f,h}p  V/eC~(P). 


The  hamiltonian  h  descends  to  h  :  P/G  — >  R  and  determines  a  reduced  dynamics  Xh  on 
P/G  by 


[/]  =  {f M P/G  VfeC°°  (P/G). 

Here  h  ([x])  =  h(x)  for  an  equivalence  class  [x]  in  P/G. 

In  what  follows,  we  work  out  the  Poisson  reduction  of  T*SE( 3)  by  SO( 3).  The  resulting 
bracket  captures  the  geometry  of  the  central  force  field  problem  studied  in  this  paper. 

Let  <  •,•  >  denote  the  pairing  between  T*SO(3 )  and  TSO( 3)  defined  by 

T*SO( 3)  x  TSO( 3)  -+  R 

<  &A  ,  WA>  =  \  tr(  oiTa  Wa  ). 


P  =  T*SE( 3)  carries  a  canonical  symplectic  structure  and  hence  a  Poisson  structure  {•,  • }p  given 

by 


{f,g}p(B,BIl,r,p) 


=  <  DBf , 


dg 

dBJI 


>  ~  <  DBg , 


df  df 
dBU  +  dr 


°l 

dp 


m  of 

dr  dp  ’ 
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where  ^  denotes  the  natural  pairing,  i.e.  the  Euclidean  inner  product  on  1R3  . 

The  group  G  =  SO(3)  acts  on  SE(3)  by  left  multiplication.  This  action  lifts  to  a  symplectic 
action  on  T*SE( 3)  given  by 

(f?,  (5,5n,r,p))  i— >  (RB,  RBTL,  Rr,  Rp). 

Thus  a  representative  for  an  equivalence  class  in  P/SO( 3)  is  given  by 

(f[  ,BTr,BTp). 

Let  A  =  BTr ,  p  =  BTp.  We  will  compute  the  reduced  Poisson  structure  on  T* S E(3) j S 0(3)  ~ 
so* (3)  X  IR3  X  IR3 .  Since  so* (3)  ~  1R3  ,  the  question  is  equivalent  to  finding  a  Poisson  structure 
on  IR3  x  IR3  x  IR3  . 

Let  f,ge  C°°( H3  x  IR3  x  H3),  and  define  f,g£  C°°(T*SE(3 ))  as 

f(B,BTl,r,p)  =  f(H,BTr,BTp). 

By  the  definition  of  reduced  Poisson  structure,  we  have 

{f,g}p/o(  n,A,/i)  =  {f,g}P(B,BU,BX,Bp). 


(The  right  hand  side  is  the  canonical  bracket  in  T*SE{3).)  Then,  by  the  canonical  bracket  on 
T*SE(  3), 


{/» 1/}P/G(n)  /O  =<  DBf,  > 

oU  11 


„  _  df  Of  dg 
<  B9:  dBU>+dr'  dp 


a fdj_ 

dr  dp ' 


Instead  of  computing  each  element  in  the  above  formula  individually,  we  compute  the  differential 
of  /.  Let  W  =  (Bdi,  B(v\fl  +  £>2), ^3,^4)  €  T^B  BB  rp^T*SE(  3).  It  generates  the  curve 

(Re**1 ,  Re*"1  (ft  +  tv2),  r  +  tv3,p  +  C  T*SE(  3) 


Thus  the  differential  is  given  by, 
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The  reduced  bracket  can  then  be  derived  as 

{/>#}p/G(n,  A,//) 


•  B 


dg_ 

dg, 


„  .  dg  -  .  dg  \  df  ^dg  ndf 

<  I  ^n-(Ax  — )-(^X  — )J  ,5— >+5—  •  5— , 

I*r  f  f[T?/ _  I<r  f nT^1 +  dJ_ 

2  1  an  an  I  2  an  an  raA  a/i  aA  a//, 


+ 


£/ 

an 


x  dg  dg 

Ax9l  +  /iX^ 


89 

an 


x  df  df ' 

AxaX  +  ,Jx^y’ 


=  _n  (<!L  x 

-l  an  anj  aA'a/i  aA '  dg 


dff  dg  dg 

+  sn'(Ax8A+,iX3f‘ 


dg 

an 


,  df  df 
Xxd\+^xdil 


df T  dfT  dfT,  (?  A  ^ 

~(an  aA  a/i  M  *  0  * 

\v  -1  0 


an  \ 

9A 

d<?  / 

d/t  / 


In  terms  of  the  notation  introduced  before,  the  matrix  form  for  the  Poisson  tensor  A  is 
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f?  :  a 

\fc  -I  0  J 

REMARK  A2.1. 

The  reduced  Poisson  structure  derived  here  is  very  closely  related  to  the  one  derived  by 
Krishnaprasad  and  Marsden  for  the  dynamics  of  a  rigid  body  with  a  flexible  attachment  [17]. 
The  key  link  is  the  geometry.  In  [17]  the  unreduced  phase  space  is  infinite  dimensional  and  is 
given  by 

P x  =  T*SO{ 3)  x  T*C 

where  C  =  {/  :  [0,  L]  — ►  1R3  |  /  is  smooth}  is  the  configuration  space  for  the  shear  beam 
attachment.  In  the  present  paper  the  unreduced  phase  space  is 

P  =  T*SE( 3)  =  T*50( 3)  x  T*I R3. 


In  both  settings,  the  reduction  is  by  50(3).  In  [17],  the  Poisson  bracket  takes  the  form 


{fy9}Poo/G  =  -  n  •  ( 


an 


df_  dg_ 
d\  d\i 


dg_  &f\ 

d\  dfi,  J 


ds 


+ 


dg_ 

d\ 


+  n  x 


dg_ 

dfj, 


ds 


x  df 

Ax  — +px 


ds, 


where  the  convected  variables  A,  /i  are  H3  valued  functions  on  [0,  L] .  As  we  let  the  flexible 
attachment  become  vanishingly  small  (L  — >  0)  with  infinitely  large  density,  {•,  •} px/o  “collapses” 
to  {' 5  ‘ }p/g  ■ 

REMARK  A2.2.  (Casimir  functions) 

A  function  /  on  a  Poisson  manifold  (P,  {•,  -}p)  is  said  to  be  a  Casimir  function  if 

{  /,  V>  }p  =  0  V  </>  6  C°°(P). 
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Note  that  a  Casimir  finction  is  automatically  a  conserved  quantity  for  any  hamiltonian  vector 
field  Xh  on  P .  If  {-,-}p  is  induced  from  a  symplectic  structure  and  P  is  connected,  then  the 
only  Casimir  functions  are  the  constant  functions.  In  the  present  context,  on  ( P/G ,  {v}p/g) 
there  are  nontrivial  Casimir  functions.  If  (II,  A,  /r  )  is  a  generic  point  on  P/G  (i.e.  a  point 
where  the  matrix  A  has  maximal  rank  (=  8)),  then  on  a  neighborhood  of  this  point  any  Casimir 
function  is  of  the  form, 

cv  =  <f>{  |  n  +  a  x  p  |2 ) 

where  (f>  :  H  — >  1R  is  an  arbitrary  real- valued  function. 

This  follows  from  the  observation  that  V  C<f,  is  in  the  null  space  of  A . 
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APPENDIX  3.  (Stability  Characterization) 

We  recall  the  Energy- Casimir  theorem  (  see  Holm,  Marsden,  Ratiu,  and  Weinstein  [16]  ). 
Consider  a  finite  dimensional  Poisson  system  with  a  hamiltonian  H 

x  =  {x,  H}(x)  =  A(x)  VH(x).  (A3.1) 

Assume  the  null  space  of  A  is  not  empty  and  is  spanned  by  VC,  where  C  is  some  (  fixed  ) 
Casimir  function.  Then  xe  is  an  equilibrium  of  (A3.1)  iff  3  A  such  that, 

V%)  =  A  VC(xe).  (A3. 2) 

Notice  that  if  C  is  a  Casimir,  so  is  any  smooth  function  C  of  C,  and  VC  also  spans  the  null 
space  of  A.  So  (A3. 2)  can  be  rewritten  as 

V(H  +  C)(xe)  =  0, 

for  a  whole  class  of  Casimirs.  Note  that  H  -j-  C  is  a  conserved  quantity  along  trajectories  of 
(A3.1).  We  have  the  following  theorem. 

THEOREM  A3.1.  (Energy- Casimir) 

If  there  exists  a  Casimir  function  C  such  that 

V(H  +  C)  (xe)  =  0,  (A3. 3a) 

(second  variation)  V2(If  +  C)(xe )  >  0,  (or  <  0)  (A3. 3b) 

then  xe  is  a  Lyapunov  stable  equilibrium  of  (A3.1). 

Proof 

Define 

V(x)  =  (H  +  C)(x)  —  (H  +  C)(xe). 
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By  assumption  V2(iJ  +  C)(xe)  is  positive  definite,  so  we  know  that  xe  is  a  strict  local  minimum. 
Thus  there  exists  a  neighborhood  U  of  xe  such  that 


V(xe)  =  0, 

V{x)  >0,  \/x£U  —  {xe}. 

Since  H  +  C  is  a  conserved  quantity  along  trajectories  of  the  given  system,  we  have  also 

V(x)  —  0,  \/x£U  —  {xe}. 


We  may  therefore  conclude  that  xe  is  Lyapunov  stable  by  a  standard  lemma  (see  e.g.  Hirsch 
&  Smale  [15],  pp.  193).  A  similar  argument  can  be  applied  in  the  case  that  V2(iT  +  C)(xe)  is 
negative  definite. 

K 

Next  we  describe  alternate  approach  to  obtain  a  stability  theorem.  Consider  the  constrained 
variational  problem 


min 

subject  to  C(x)  =  b, 


(A3. 4) 


where  b  is  a  constant  representing  prescribed  data.  The  Lagrangian  corresponding  to  this 
optimization  problem  can  be  written  as 

L(x,X)  =  H{x)-\C(x),  (A3. 5) 


with  A  G  1R.  The  first  order  necessary  conditions  for  (A3. 4)  then  coincide  with  (A3. 2).  We  now 
recall  the  following  lemma. 

LEMMA.  (  See  e.g.  Bertsekas  [4],  p.  68  ) 

Let  P  be  a  symmetric  matrix  and  Q  a  positive  semidefinite  symmetric  matrix,  both  of 
dimension  n  X  n.  Assume  that 

(a;,  Px)  >0,  V  x  0,  with  {x,  Qx)  =  0, 
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then  there  exists  a  (large,  positive)  scalar  a  such  that 

P  +  oi  Q  >  0, 

i.e.  P  +  aQ  is  positive  definite. 

I 

We  can  now  state  the  stability  criterion  as  follows. 

THEOREM  A3. 2. 

Suppose  that  xe  and  Ae  €  1R  are  such  that 

\7xL(xe,Xe)  =  0,  (A3. 6a) 


and,  moreover, 

(h,  V^.L(xe,  Xe)  h)  >  0,  V  h  0  with  (VC'(xe),  h)  =  0.  (A3. 66) 


Then  xe  is  a  Lyapunov  stable  equilibrium  of  (A3.1). 

Proof 

Let 

P  —  V^,L(.re,  Ae), 

Q  =  VC(xc)  VC(re)T, 

so  that  by  hypothesis  P  and  Q  satisfy  the  conditions  of  the  previous  lemma.  Thus  we  can  find 
a  E  IR  such  that  P  +  aQ  is  a  positive  definite  matrix.  Now,  with  the  notation  b  =  C(xe),  define 
the  augmented  Lagrangian  by, 

La(x,  A)  =  H(x)  —  XC(x)  +  —a  ^C(x)  —  bj  . 


Then 
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VxLQ(xe,  Ae)  =  VH(xe)  +  \eVC(xe)  +  a  (c(xe)  -  b)VC(xe)  = 

V2xLa(xe,  Ae)  =  V2Jff(xe)+  (Ae  +  n  (d(xe)-bj)V2C(xe)  +  a  VC(xe)VC(xe)T 

—  P  +  a  Q  >  0. 

Thus  the  augmented  Lagrangian  satisfies  the  requirement  of  a  Lyapunov  function  and  Theo¬ 
rem  A3.1  can  be  applied  to  conclude  that  xe  is  a  Lyapunov  stable  equilibrium  of  system  (A3.1). 

I 

REMARKS. 

(a)  Conditions  (A3. 6)  form  a  set  of  sufficient  conditions  for  xe  to  be  a  constrained  local 
minimizer  of  (A3. 4). 

(b)  In  application  of  Theorem  A3.1,  we  would  search  for  a  suitable  Casimir  C  to  fulfill  the 
condition  (A3. 3).  However,  in  application  of  Theorem  A3. 2  we  can  fix  a  particular  Casimir 
C  and  a  scalar  A  satisfying  (A3. 6a),  and  then  attempt  to  verify  (A3. 6b).  The  analysis  in 
Sections  7.2  and  7.3  illustrate  the  differences  between  the  two  schemes. 

(c)  With  appropriate  hypotheses,  both  Theorem  A3.1  and  A3. 2  can  be  generalized  to  cases  in 
which  there  are  n  independent  Casimirs,  and  the  underlying  space  is  infinite  dimensional. 
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